A charged particle in a magnetic field - Jarzynski Equality 
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Abstract: We describe some solvable models which illustrate the Jarzynski theorem and related 
fluctuation theorems. We consider a charged particle in the presence of magnetic field in a two 
dimensional harmonic well. In the first case the centre of the harmonic potential is translated with 
a uniform velocity, while in the other case the particle is subjected to an ac force. We show that 
Jarzynski identity complements Bohr-van Leeuwen theorem on the absence of diamagnetism in 
equilibrium classical system. 

PACS numbers: 05.70.Ln, 05.40. a, 05.40.Jc 
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Most processes that occur in nature are far from 
equilibrium and hence cannot be treated within the 
framework of classical thermodynamics. The traditional 
nonequilibrium statistical mechanics deals with systems 
near equilibrium in the linear response regime. Its suc- 
cess has lead to the formulation of fluctuation-dissipation 
relation. Onsagar's reciprocity relations and the Kubo- 
Grecn formulae, etc. However, very recent develop- 
ments in nonequilibrium statistical mechanics have re- 
sulted in the discovery of some exact theoretical re- 
sults for systems driven far away from equilibrium and 
are collectively called fluctuation theorems[l]. These re- 
sults include entropy production theorems [2], Jarzyn- 
ski equality [3], Crooks relations[4], Hatano-Sasa iden- 
tity [5], etc. Some of the above relations have been verified 
experimentally on single nanosize systems in physical en- 
virnoments where fluctuations play a dominant role [6, 7]. 

The concept of free energy is of central importance 
in statistical mechanics and thermodynamics. With the 
help of free energy one can calculate all the phases of 
a system and their physical properties. However, the 
free energy of the system relative to an arbitrary refer- 
ence state is often difficult to determine. Jarzynski equal- 
ity (JE) relates non-equilibrium quantities with equilib- 
rium free energies [3] . In this prescription, initially the 
system is assumed to be in equilibrium state determined 
by a thermodynamic parameter A defined by a control 
parameter A,4 and is kept in contact with a heat bath at 
temperature T. The nonequilibrium process is obtained 
by changing the thermodynamic control parameter A in a 
finite time r according to a prescribed protocol A(i), from 
\a = A(t=0) to some final value A^ = X(t = t). The 
final state of the system at time r (at the end of the pro- 
tocol) will in general not be at equilibrium. It will equi- 
librate to a final state B(= Xb) if it is further allowed to 
evolve by keeping parameter Xb fixed. JE states that 
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Where AF is the free energy difference between equi- 
librium states A and B. Angular bracket (...) denotes 
the average taken over different realizations for fixed pro- 



tocol A(t). W is work expended during each repeta- 
tion of the protocol and is a realization dependent ran- 
dom variable. Jarzynski's theorem has been derived us- 
ing various methods with different system dynamics [3- 
5,8]. This remarkable identity provides a practical tool 
to determine equilibrium thermodynamic potentials from 
processes carried out arbitrarily far away from equilib- 
rium. This identity has been used to extract equilibrium 
free energy differences in experiments. Work distributions 
have been calculated analytically for several model sys- 
tems and tested against various fluctuation theorems [9- 
11]. JE has been generalized to arbitrary transitions 
between nonequilibrium steady states by Hatano and 
Sasa[5], which has also been verified experimentally [7]. 

In our present work, we determine the work distri- 
butions analytically for two different models. In both 
these cases charged particle dynamics in a two dimen- 
sional harmonic trap in the presence of a magnetic field 
is considered. In the first case (i), centre of the harmonic 
trap is dragged with a uniform velocity whereas in the 
case (ii), the particle is subjected to an ac force. We 
show that JE is consistent with Bohr-van Leeuwen theo- 
rem. We also discuss steady state fluctuation theorems 
[SSFT] and energy loss in driven systems. 

We consider a charged particle motion in two dimen- 
sions (x-y plane) in the presence of a time dependant 
potential U(= U(x,y,t)). An external magnetic field 
B along z-direction will produce Lorentz force on the 
charged particle. The interaction of the particle with en- 
vironment can be treated via frictional force along with 
concomitant fluctuations. The appropriate equations of 
motion are given by the Langevin's equations [12, 13] 
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where the random force field is a Gaussian white 

noise. i.c, 



(Ut)&(t ))=D8 a0 5(t-t), 
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with a, (3 = x,y, e is the charge of the particle. D = 
2jksT is a consistency condition for the system to 
approach equilibrium in the absence of time independent 
potential . The friction coefficient is denoted by 7. 

This problem (for time independent potentials) was 
earlier considered [13] to elucidate the crucial role played 
by the boundary conditions in the celebrated theorem 
of Bohr-van Leeuwen on the absence of diamagnetism 
in classical systems [13, 14]. This theorem states that in 
equilibrium for classical systems the free energy is inde- 
pendent of magnetic field. Hence, diamagnetism does 
not exist in classical statistical mechanics. 

We restrict our analysis to overdamped regime where 
the corresponding dynamical equations become 
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The associated Fokker-Planck equation [15] leads to an 

_(U1e±e1) 

equilibrium distribution P e oc e k B T in a stationary 
regime for time independent potential U(x,y). This dis- 
tribution P e is independent of magnetic field consistent 
with Bohr-van Leeuwen theorem. 

We consider in this work two cases for time dependent 
potentials. For case(i) U(x,y,t) = \k\r — r* (t)\ 2 , where 
f is a two dimensional vector r = ix + jy and r* (t) = 
vt(i i and j are unit vectors along x and y direction 
respectively. Here the particle is in a harmonic poten- 
tial whose centre is dragged along with an uniform speed 
in a diagonal direction. This problem can also be 
solved for the motion of the centre in an arbitrary direc- 
tion with a different protocol. For case(ii), U(x,y,t) = 
\k[x 2 + y 2 ) — Ax sin(tJt) . Here the particle in a two di- 
mensional harmonic well is subjected to an ac force in 
x-direction. This problem can also be solved for ac driv- 
ings in both x and y direction with different amplitudes 
and with a phase difference. 

We rewrite equations (5) and (6) using the variable[12] 
z = x + iy,(i = ^/—l). For the case(i) we get 
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where p = £(t) =&(*)+ if „(*)> <?*(*) = vt(l + i) 

and C = £M 

For the case(ii), we get 



-kpz p 



7 



7 



(Asin(wi)+f(i)). 



(8) 



Thermodynamic work W done on the system by an ex- 
ternal agent during a time interval t is given by [9 ,10], for 



casc(i), 



W = -kv {(x - vt) + (y - vt )}dt , (9) 



and for case(ii) 



W = -Alu / cos(w< )x(t )dt, 
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The formal solution of equations (7) and (8) respectively 
are 



k p 



k 



z(t) = ^ocxp( — pt) H — / dt exp( — pit — t )){kg*(t ) 
7 7 Jo 7 



and 
z(t) 



<(t )}, 



z exp( pt) 

7 

+Asm(ujt))}. 



(11) 



7 Jo 



t , k , 

dt exp( — p(t - t )){£(< ) 



7 



(12) 



where zq = xq + iyoi %o and yo are initial co- 
ordinates of the particle. The initial distribution for Xq 
and yo is assumed to be equilibrium canonical distribu- 



It may be readily 



tion P e (x , yo,*) = cx p[" 2 
noted from equations (10), (11), (12) and (13) that work 
done as well as particle co-ordinates at later times are lin- 
ear functional of Gaussian variables £ x (t) and £ y (t) and 
hence their distributions are Gaussian. We calculate full 
probability distribution for W for both cases following 
closely the procedures adpoted in refs[9,10] . Without go- 
ing into further details of algebra we give our final results 
which will be further analyzed. 

For the case(i), the average work done (W) is 

(W) = 2~/v 2 {t - ^-(1 - exp(-fc*i) cos(Qt)) - — - sin (fit) 
k k 



cxp(-/c*<)} - 7v 2 2C{-sin(fit)exp(-fc*i) 
k 

(1 - exp(-Fi) cos(fit))}, 
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where fi = ^r^^rj an d k* = -7^^ ■ The above equa- 
tion(13) agrees with the result obtained in refs[9,10] for 
B = 0. The variance of the work is given by 
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where f3 = -j^f- the full probability distribution P(W) 



is 



P(W) = 
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JE given in equation(l) follows immediately from the 
above expression, with 
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The equation(16) implies AF = 0, indicating that the 
equilibrium free energy of a particle in a harmonic po- 
tential is independent of magnetic field, consistent with 
the Bohr-van Leeuwen theorem. Needless to say that in 
the present case the free energy of the oscillator is also 
independent of position of the centre of the harmonic po- 
tential as expected on general grounds. However, it is in- 
teresting to note that the thermodynamic work W in the 
transient state depends explicitly on the magnetic field 
B. In the presence of magnetic field relaxation rate of the 
system r r (= 7 ^ c ^ ) depends on the magnetic field. It 
increases with the strength of the magnetic field. Hence 
the magnetic field gives an additional control over the 
relaxation time in experimental situation to verify the 
above and subsequent results. 

Only in the asymptotic time limit t -> oo(i » 
r r ), where the system docs not retain the memory of 
the initial state, one obtains for the averaged work done 
(W s ) in a system in time interval t, 
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Hence power (p) delivered to the system in the steady 
state is constant (p = 2jv 2 ). In this state thermodynamic 
work is essentially a mechanical work delivered to the sys- 
tem by a moving trap with a speed s/2v along the diago- 
nal direction in x-y plane. The particle in this state set- 
tles to a Gaussian distribution with the same dispersion 
as in the case of canonical equilibrium distribution. How- 
ever, the centre of the position of the particle distribution 
lags behind the instantaneous minimum of the confin- 
ing potential by a distance I = u 1 along the diago- 
nal line. The harmonic potential pulls the particle with 
a force kl, at speed \[2v. Thus the power delivered is 
kly/2v = 2jv 2 . This power is dissipated as heat into the 
surrounding medium. In this regime, the fluctuations in 
work (W s ) is 



(w 2 ) = (w s ) 2 + \(w s 
p 
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In the steady state distribution of work W s obeys the 
SSFT[9] namely 



P(W S ) 



P(-w.) - < w > 

It is important to note that both work and its distri- 
bution in the steady state do not depend on the mag- 
netic field. Hence the experimentally obtained results[7] 
to verify Hatano-Sasa identity [5] for transition between 
noncquilibrium steady states, remains unaltered irre- 
spective of the magnetic field being present or not . Here 
transition between noncquilibrium states are induced by 
varying the speed of the trap [7]. 

We now turn to case(ii). As mentioned earlier for 
this case the work distribution is Gaussian and is char- 
acterized completely by the first moment and the vari- 
ance. Expression for (W) is too lengthy and is given in 



appendix A. The variance is given by 



(W 2 ) - (W) 2 = -((W)-AF) 



(20) 



Where AF = sin 2 (car). Using the distribution of 
work one can readily verify the JE, namely (e~@ w ) = 
e -/3AF^ - 1S foe free energy difference between two 
thermodynamic states. The initial (t = 0) and final 
states (t = t) are characterized by a two dimensional 
harmonic potential with an additional tilt of magnitude 
zero and — Axsin(uiT), respectively. Note that AF is in- 
dependent of magnetic field, again reassuring Bohr-van 
Leeuwen theorem. However, the work distribution is an 
explicit function of magnetic field. 

We now concentrate on the statistics of the work done 
W s in the asymptotic regime. In this regime probabil- 
ity distributions are time periodic with a period 2j. We 
calculate averaged work done over one period , which 
is 
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Similarly variance (V s ) of the work averaged over a 
period of oscillation in is given by 



(V s ) = (W 2 ) - (W 8 f 



(Vs) = -(W s ) 



(23) 
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In the time periodic state averaged input energy is dis- 
sipated into system as a heat [9]. Thus one can identify 
(W s ) as a hysteresis loss in the medium. Since the prob- 
lem being linear we find that the time averaged hystere- 
sis loss is independent of temperature. However, it de- 
pends explicitly on the magnetic field and is a symmet- 
ric function of magnetic field. Thus the magnetic field 
becomes a relevant variable in the time periodic asymp- 
totic state. However, it must be noted that variance in 
the input energy can not be identified with the heat fluc- 
tuations [9]. In this time periodic state work done W s over 
a period satisfies the SSFT, namely, equation(19). 

It may be noted that the validity of SSFT for work 
done over a single period (equation(24)) is restricted 
only to overdamped linear models as in the present 
case. In general this will not hold good in nonlinear situa- 
tions[17] . However, one can show that if one instead con- 
siders work done over a large number of periods, indeed 
SSFT holds even for nonlinear models. The convergence 
of SSFT on accessible time scales have been studied in the 
previous literature[18]. We have calculated separately 
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the work distribution for W s in the inertial regime which 
satisfies SSFT. Moreover we have also shown that one 
can obtain the orbital magnetic moment in this nonequi- 
librium state without violating the Bohr-van Leeuwen 
theorem. These results will be published elsewhere[16]. 

In conclusion we have solved analytically the work dis- 
tribution of a charged particle in the presence of magnetic 
field in two different cases. The first being the case where 
minimum of the harmonic potential is dragged with an 
uniform velocity and in the second case particle is sub- 
jected to an ac force. For both cases JE is verified and this 
equality complements Bohr-van Leeuwen theorem on the 
absence of diamagnetism in classical system. In case(i) we 
have shown that such a distribution in steady state does 



not depend on magnetic field and satisfies SSFT. As op- 
posed to this case in a time periodic assymptotic state for 
casc(ii) magnetic field becomes a relevant variable. The 
hysteresis loss over a cycle depends explicitly on the mag- 
netic field. Relaxation time in our system can be con- 
trolled by magnetic field. All our results are amenable to 
experimental verification with charged beads in magnetic 
field. 



APPENDIX A: 

An expression for (W) 
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where 
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oi = -2A 2 uj 2 (k*' +n 2 +lu 2 ), 

a 2 = k*{k 2 + lo 2 (1 + C 2 )} + nC(k 2 + iio 2 - lo 2 C 2 ), 
ci = (k* 2 + lu 2 (1 - C 2 )) 2 + AC 2 u\ 

C 2 = 2[k* 2 + (0 - L0) 2 }[k* 2 +(n + Uj] 2 }. 
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In the absence of magnetic field, (W) reduces to 



(W) = 



A 2 k 



2 7 (fc 2 + cj 2 ) 
A 2 u 2 t 



sin 2 (wt) 



A 2 uj 



47 (uj 2 + k 2 ) 
A 2 uj 2 k A 2 uj 2 



sin(2w<:) 



2 1 {k 2 +uj 2 ) 1 (k 2 +uj 2 ) 2 7 (fc 2 + w 2 ) 2 
exp(— kt)[k cos(ujt) — lu sin(utf)] (27) 
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